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[i]  Euphotic  zone  depth,  zl%,  reflects  the  depth  where  photosynthetic  available  radiation 
(PAR)  is  1%  of  its  surface  value.  The  value  of  z1%  is  a  measure  of  water  clarity,  which  is 
an  important  parameter  regarding  ecosystems.  Based  on  the  Case-1  water  assumption, 
z1%  can  be  estimated  empirically  from  the  remotely  derived  concentration  of  chlorophyll-a 
([Chi]),  commonly  retrieved  by  employing  band  ratios  of  remote  sensing  reflectance  (R^). 

Recently,  a  model  based  on  water’s  inherent  optical  properties  (IOPs)  has  been  developed 
to  describe  the  vertical  attenuation  of  visible  solar  radiation.  Since  IOPs  can  be  near- 
analytically  calculated  from  so  too  can  Zj%.  In  this  study,  for  measurements  made  over 
three  different  regions  and  at  different  seasons  (z]%  were  in  a  range  of  4.3-82.0  m  with 
[Chi]  ranging  from  0.07  to  49.4  mg/m3),  z1D/o  calculated  from  was  compared  with 
z1%  from  in  situ  measured  PAR  profiles.  It  is  found  that  the  zl%  values  calculated  via 
/fns-derived  IOPs  are,  on  average,  within  ~14%  of  the  measured  values,  and  similar  results 
were  obtained  for  depths  of  10%  and  50%  of  surface  PAR.  In  comparison,  however,  the 
error  was  ~33%  when  ziyB  is  calculated  via  tffS-derived  [Chi].  Further,  the  importance  of 
deriving  euphotic  zone  depth  from  satellite  ocean-color  remote  sensing  is  discussed. 
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L  Introduction 

[z]  Euphotic  zone  depth,  Zi%,  reflects  the  depth  where 
only  1%  of  the  surface  photosynthetic  available  radiation 
(PAR)  remains  [A7rA,  1994].  zl%  is  a  measure  of  water 
clarity,  which  is  not  only  a  quality  index  of  an  ecosystem 
but  also  an  important  property  for  primary  production 
[Rehrenfeld  and  Falkowski,  1 997;  Platt  and  Sathyendranath, 
1988]  and  heat  transfer  [Chang  and  Dickey,  2004;  Kara  et  al, 
2005;  Sathyendranath  et  al,  1991]  in  the  upper  water 
column.  Based  on  Case-1  water  assumptions  [Morel, 
1988;  Morel  and  Prieur T  1977],  zl%  can  be  estimated  from 
remotely  derived  concentration  of  chlorophyil-a  or  total 
pigments,  with  either  a  spectrally  integrated  form  [Morel, 
1988]  or  a  spectrally  resolved  formalism  [Prieur  and 
Sathyendranath,  1981;  Sathyendranath  et  al,  1989a]. 

[j]  Case-1  waters  are  those  whose  inherent  optical  prop¬ 
erties  [Preisendorfer,  1976]  can  be  adequately  described  by 
phytoplankton  (represented  by  chlorophyll  concentration) 
[Gordon  and  Morel ,  1983;  IOCCG,  2000;  Morel,  1988; 
Morel  and  Prieur ,  1977].  Case-1  waters  thus  require  the 
optical  properties  of  other  optically  active  constituents  (such 
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as  Colored  Dissolved  Organic  Matter  and  suspended  par¬ 
ticles)  closely  follow  the  optical  properties  of  phytoplankton 
[Morel,  1988;  Morel  and  Maritorena,  2001].  Because  the 
definition  of  Case- 1  water  is  not  based  on  the  geographical 
location,  nor  based  on  the  value  of  chlorophyll,  it  is  difficult 
to  know  a  priori  if  a  water  body  fits  the  Case- 1  definition 
when  it  is  measured  by  a  remote  sensor.  Frequently  open 
ocean  waters  are  assumed  as  Case-1,  but  recent  studies 
[Lee  and  Hu,  2006;  Mobley  et  al ,  2004]  have  shown  that 
tliis  could  be  problematic,  as  in  general  water’s  optical 
properties  are  not  determined  by  phytoplankton  alone 
[Sathyendranath  et  al,  1989b]. 

[4]  To  incorporate  the  dependence  of  subsurface  light 
field  on  other  water  constituents  (such  as  CDOM),  Liu  et  ai 
[2002]  developed  a  numerical  model  (via  look-up-table) 
that  uses  information  about  chlorophyll  concentration, 
CDOM  absorption  and  particle  scattering  coefficient  as 
inputs  to  describe  the  vertical  distribution  of  dowuwelling 
irradiance.  This  approach,  similar  as  those  based  on  Case- 1 
assumption  [Morel  and  Antoine,  1 994;  Ohlmann  and  Siegel, 
2000;  Sathyendranath  et  al,  1989a],  however,  still  requires 
accurate  information  about  chlorophyll  concentration  when 
applied  to  ocean-color  remote  sensing.  Liu  et  al  [2006]  are 
successful  using  this  approach  for  waters  of  the  West 
Florida  Shelf, 

[5]  In  another  empirical  approach,  Mueller  and  Lange 
[1989]  developed  relationships  to  calculate  zy%  directly 
from  the  diffuse  attenuation  coefficient  at  490  nm.  The 
associated  empirical  constants,  however,  were  developed 
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for  a  specific  region  and  specific  temporal  ranges  [Mueller 
and  Lange ,  1989],  it  is  then  difficult  to  apply  the  empirical 
relationships  to  the  global  oceans  at  different  seasons. 

[^]  Rccenlly,  based  on  the  principle  that  the  vertical 
variation  of  subsurface  light  field  is  determined  by  inherent 
optical  properties  (IGPs),  Lee  et  at  [2005b]  developed  an 
analytical  model  to  describe  the  vertical  attenuation  of 
down  welling  vector  irradiance  in  Ihe  visible  domain  (£Vis, 
350-700  nm).  In  this  IOP-centercd  approach,  the  vertical 
attenuation  coefficient  of  £v IS  (J^vis(^))  is  modeled  as  an 
analytical  function  of  waler's  absorption  and  backscattering 
coefficients.  For  vertical  distributions  of  Ey]S(z)  ranging 
from  0. 1  to  1 00%  of  the  surface  value,  the  modeled  £v  IS(z) 
is  accurate  to  within  ^6%  of  the  actual  value  [Lee  et  al , 
2005b],  for  a  data  set  simulated  using  Hydrolight  [Mobley, 
1995].  Since  the  vertical  distribution  of  £Vis(z)  is  nearly 
identical  to  that  of  B4/2(z)  [Lee  et  al *  2005  b;  Morel  and 
Gentili,  2004],  zJ%  can  be  easily  calculated  when  lOPs  (the 
absorption  and  hackscattering  coefficients  at  490  nm,  in 
particular)  are  known,  either  from  in  situ  measurements  or 
from  remote  sensing  of  ocean  color. 

[7]  Regardless  of  whether  one  uses  the  Case- 1 
approaches  or  the  recent  lOP-cenlered  approach,  there  has 
been  no  test  or  validation  of  z\%  algorithms  with  a  broad- 
range  in  situ  data.  On  the  other  hand,  global  or  basin  scale 
estimation  of  from  ocean  color  has  been  incorporated 
into  studies  of  new  production  [Behrenfeld  et  al ,  2005]  and 
the  Sverdrup  hypothesis  [Siegel  et  al,  2002].  To  ensure 
reliable  products  for  these  biogeo  chemical  studies,  it  is 
desirable  to  know  the  accuracy  of  remotely  derived  zl% 
for  broad  range  of  waters. 

[s]  In  this  study,  for  measurements  made  in  the  Arabian 
Sea,  the  Monterey  Bay  and  the  Gulf  of  Mexico  that  covered 
both  oceanic  and  coastal  waters  and  measured  at  different 
seasons,  zj%  (and  z10%1  are  calculated  from  spectral 
remote-sensing  reflectance  using  both  the  Case- 1  and  IOP- 
centered  approaches.  With  an  emphasis  on  computational 
efficiency,  only  spectrally  integrated  approaches  arc  studied 
here.  The  calculated  zT%  (jr  represents  I,  10,  or  50  in  this 
article)  values  are  then  compared  with  those  from  profiles  of 
PAR  measurements  to  evaluate  the  performance  of  the 
approaches.  Further,  we  discussed  the  unique  characteristics 
of  property  ri%  (and/or  Zio%)  and  its  importance  in  ocean 
color  remote  sensing, 

2*  Data  and  Methods 

[9]  Field  measurements  from  three  different  regions  are 
used  for  testing  and  evaluating  the  methods  for  deriving  zxVt 
from  Rn.  These  regions  are:  the  Monterey  Bay  (data 
collected  in  September  and  October  19S9),  the  Gulf  of 
Mexico  (April  and  June  1993)*  and  the  Arabian  Sea 
(December  1994).  These  water  environments  include 
open-ocean  oligotrophic  waters  (the  Loop  Current,  and 
Arabian  Sea),  coastal  high-productivity  waters  (the 
Monterey  Bay,  the  West  Florida  Shelf),  and  turbid  Mis¬ 
sissippi  River  plume  waters.  The  range  of  chlorophylls 
concentration  measured  was  ^0, 07-49.4  mg/m*  with  zi% 
ranging  from  ^4.3  to  82.0  m.  During  these  measurements 
the  above-surface  solar  zenith  angles  (0n)  were  between 
8°  and  80°  (30%  of  them  with  0a  >  50°).  For  the  conditions 


with  the  Sun  covered  by  clouds,  is  taken  in  our  calcu¬ 
lations  as  45°  as  in  Sathyendranath  et  a l  [1989a]. 

[10]  There  are  65  stations  that  include  measurements  of 
both  remote-sensing  reflectance  and  vertical  profiles  of  PAR. 
Remote-sensing  reflectance  at  wavelength  A,  A),  was 
calculated  from  measurements  made  above  the  sea  surface  as 
described  in  Carder  and  Steward  [1985],  with  upwelling 
radiance,  downwelling  irradiance,  and  down  welling  sky 
radiance  measured  by  a  handheld  spectroradiometer.  Tire 
methodology  of  determining  in  situ  Rr,{ A)  is  described  in 
detail  in  NASA  protocols  [Mueller  ct  al,  2002a]. 

[11]  Instantaneous  PAR  (400-700  nm)  in  the  upper  water 
column  at  time  t  {PARt{z))  was  measured  by  lowering  a  PAR 
sensor  (Biospherical  Instruments,  Inc.)  from  surface  to  depth 
(z  m,  positive  downward),  z  was  measured  by  a  pressure 
sensor,  and  was  accurate  to  within  0.  t  m  with  a  fall  rate  of  1  m 
per  second.  During  the  deployment  of  the  PAR  sensor  in  the 
water  column,  a  deck  cell  was  operated  10  measure  simulta¬ 
neously  the  above-surface  downwelling  irradiance  at  time 
t  (£j(490,  /)),  in  order  to  correct  for  any  variations  of  input 
solar  light  resulting  from  passing  clouds.  The  vertical 
profile  of  PAR  free  of  cloud  effects  is  [Smith  et  al,  1984] 


PAR(z)  =  PARt{z) 


£,(490  nm,  0) 
£5(490  nm,  0  + 


0} 


From  these  PAR(z)  values,  the  ratio  of  PAR  at  depth  to 
surface  PAR  (/J47?(0))  was  calculated  as 


,  ^  PAR[z) 

rpAR{z)  -  pAR{QY  (2) 

^jnv.t  ’10%  and  z,%  (corresponding  to  rMK  equals  0.5,  0.1 
and  0,01,  respectively)  were  then  determined  from  the 
vertical  profile  of  rPAR{z),  Tt  is  rare  to  have  recorded  depth 
with  rPAR (z)  exactly  die  desired  ratio  (10,0%,  for  example)* 
The  measured  zv%  value  is  then  an  approximation  by 
exponentially  interpolating  rPAR(z)  between  ~LLc%  and 
0.9jt%.  Because  the  depth  interval  reported  for  the  vertical 
profiler  is  less  than  1  m,  the  maximum  error  in  zAyt  is  less 
than  0.5  m. 

[n]  For  some  earlier  measurements  (the  9  Monterey  Bay 
stations  in  1989  and  10  stations  in  the  Gulf  of  Mexico  in 
1 993),  the  PAR  sensor  was  not  sensitive  enough  to  provide 
a  reading  when  PAR  is  1%  of  PAR(Q ).  For  those  stations, 
only  zI0%  and  z5^  were  directly  obtained  from  rPAR(z) 
profiles. 

[n]  For  the  other  stations  where  both  zi0%  and  z!%  could 
be  obtained  from  rPAR[z),  it  is  found  that  a  good  linear 
relationship  exists  between  and  rl%  (see  Figure  1 ).  For 
Z]%  in  the  range  of  4.3-64.4  m  (zi0%  in  the  range  of  2. 1  - 
29.0  m)t  by  forcing  the  linear-regression  line  to  pass 
through  the  origin,  the  following  relation  was  obtained, 


z1%  =2.25z|0%,  (3) 

with  a  coefficient  of  determination  of  0.99  (n  =  43).  In  this 
regression  analysis,  the  two  measurements  with  of 
^82.0  m  (the  black  points  in  Figure  1)  were  excluded, 
simply  because  they  were  significantly  beyond  the  general 
linear  trend  with  too  few  points  in  that  range  to  make  a 
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Figure  h  Relationship  between  zJ%  and  zt0%  from  in  situ  PAR  profiles.  The  two  black  dots  were 
excluded  in  the  linear  regression,  and  the  range  of  zl0%  is  2 A  -29.0  tn. 


confident  generalization.  For  the  same  ranges  of  zl%  and 
zh,%>  simulations  using  Hydrolight  [Mobky>  1995]  with 
different  sun  angles  and  lOPs  indicate  a  nearly  identical 
relationship  as  equation  (3)  (with  a  slope  of  2.21).  These 
results  then  suggest  that  equation  (3)  does  represent  a 
reliable  relationship  between  zl%  and  zj0%,  at  least  for  zJ0% 
in  this  specific  range  of  —2  —  30  m.  Therefore,  for  those 
stations  where  no  =\ym  were  obtainable  from  rpAH(z)  and  zw% 
was  in  the  range  of  2-30  m,  zy%  values  representing  PAR 
profiles  were  derived  by  applying  equation  (3).  tn  essence, 
in  this  study,  equation  (3)  is  only  utilized  for  range- 
constrained  interpolation,  not  for  extrapolation, 

[u]  Usually,  zw%  is  assumed  to  be  half  of  zr/a  (zw%  is 
called  midpoint  depth  in  Kirk  [1994]).  This  is  based  on  the 
assumption  that  the  attenuation  coefficient  for  PAR  (^par) 
does  not  change  with  depth.  However,  because  PAR  con¬ 
verges  to  wavelengths  with  less  attenuation  coefficients 
when  PAR  propagates  from  surface  to  deeper  depths,  /fPAR 
is  always  smaller  at  depth  than  at  surface  (the  difference  can 
be  a  factor  of  2  or  3)  even  for  a  vertically  homogeneous 
water  environment  [Lee  et  ai ,  2005b;  More /,  1988;  Zaneveld 
et  oLt  1993],  and  thus  zWyu  is  shallower  than  half  of  zr/a. 

3,  Derive  zi%  From  /frs 
3,L  IQP-Centeml  Approach 
[is]  The  vertical  distribution  of  down  welling  trrudtance 
in  the  visible  domain  (£VIS,  350-700  run)  is  described  as 

Em(z)  =  (4) 

Kvrsfc),  the  attenuation  coefficient  of  £Vis»  however,  is  no 
longer  treated  as  independent  of  z  but  a  function  of  both  z 
and  inherent  optical  properties  [Lee  el  ai ,  2005b], 

*«.(•) (5) 


with 

|  /M  -  [,Vo  +  Vi  (o(4£>0)  )o  s+  x-3^(490)](I  +  a0sin(^)), 

l  ^2  =  [Co  +  (|fl(490)  +  (ih(«0)]{oi  +  cuK)), 

where  #(490)  and  &h(49Q)  are  water  absorption  and 
backscattering  coefficients  at  490  nra,  and  Xoa.2>  Co,u  and 
V(Qjt2  are  model  constants  derived  from  Hydrolighl  simula¬ 
tions  [Lee  el  ai ,  2005b].  Given  a  sun  angle  and  values  of 
a(490)  and  6b(490),  Avjsfc)  can  then  be  calculated  for  any 
depth. 

[ie>]  Defining  the  optical  depth,  T£,  for  £Vis  as 


Kns[?)z  ~  r£f 


(?) 


then  the  dcplh  where  £Vis{z)  is  l%of£Vis(0)  is  -ln(0,01)  = 
Tg  =  4.605.  Combining  equations  (5)  and  (7),  after  simple 
math  manipulations,  a  cubic-polynomial  equation  with  z  as 
the  variable  is  obtained 

z3  Tyjz2  +.Viz+yi  =  0,  (8) 


with  3  functions  of  Ku  K2  and  rEt 


y  i 


>2 


m 


Mathematically  there  are  three  solutions  (one  negative,  two 
positives)  that  satisfy  equation  (8),  but  it  is  the  smaller, 
positive,  value  that  is  consistent  with  radiative  transfer 
theory.  In  this  study,  as  examples  to  demonstrate  this  IOP- 
centered  approach,  z  for  rE  equals  4.605  (£Vist)^vis(0)  = 
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1%),  2.303  (10%),  and  0.693  (50%)  are  calculated  for  those 
measurements,  respectively. 

[i?]  Before  the  derivation  of  z  for  different  r&  inherent 
optical  properties  (ri(490)  and  6b(490)  in  particular)  were 
derived  from  /?rs  using  the  quasi-ana  lyrical  algorithm 
(QAA)  [Lee  et  ai *  2002].  QAA  is  an  algorithm  for  deriving 
absorption  and  bae Scattering  coefficients  of  optically  deep 
waters  from  occan-color  remote  sensing.  Its  concept  and 
architecture  are  documented  in  detail  in  Lee  et  ai  [2002]. 
Here  a  slightly  updated  version  of  QAA  (v4)  was  applied, 
and  the  details  of  the  updates  are  provided  in  Appendix  A. 

[tx]  It  is  necessary  io  point  out  that  EVis  m  equation  (4) 
represents  downwdling  vector  iiradiance  in  the  range  of 
350-700  nm,  and  is  measured  by  energy  (W/m2),  whereas 
the  PAR  sensor  measures  both  downwclling  and  upwelling 
scalar  irradiance  in  the  range  of  400-700  nm  and  ts 
measured  by  the  amount  of  photons  (quanta/m2/s).  For  die 
same  wavelength  range  (400-700  nm),  Hydro! jght  simu¬ 
lations  [Morel  and  Gentili,  2004]  indicate  that  z1%  is  about 
the  same  when  it  is  measured  either  by  £Vis  or  by  PAR 
(Table  2  of  Morel  and  Genii l i  [2004]).  Also,  sensitivity  tests 
with  Hydro! ight  simulations  indicate  that  the  vertical  profile 
of  £vis  is  nearly  identical  to  that  of  PAR,  partially  because 
upwelling  irradiance  is  in  the  order  of  ^5%  of  downwdling 
irradiance,  and  that  the  irradiance  in  the  range  of  350— 
400  nm  is  small.  In  general,  the  difference  between  zx% 
measured  by  either  EVi s  or  PAR  is  small  (less  than  10%). 
Therefore,  z*%  calculated  from  equation  (4)  is  considered 
equivalent  and  comparable  to  that  determined  from  the 
vertical  profiles  of  /34tf(z). 

3.2.  |€hl]-Centerecl  Approach 

[id]  Based  on  the  Case- 1  water  assumption,  and  in  a 
spectrally  integrated  form,  zj%  can  be  estimated  from 
chlorophyll -a  concentration  ([Chi])  of  the  surface  layer 
(A.  Morel,  personal  communication,  after  the  statistical 
analysis  of  the  shape  of  [Chi]  vertical  profiles  [Morel  and 
Berthon,  1989]), 


*,*=34.0  ([CA/])-*”  (10) 


Because  surface  [Chi]  can  be  calculated  from  ^(A)  by  the 
operational  OC4v4  algorithm  [O  'Reilly  et  ai ,  2000], 


[Chi\  =  10 


0  366— 3  067/74  I  np 4 0  MQ/?-  I  532/ 


00 


with 


log( 


max(/Z„(440, 490,510}) 
R.A5 55) 


)■ 


(12) 


it  is  then  quite  straightforward  to  estimate  %i%  via  tills  route. 


4.  Results  and  Discussion 

[2a]  To  evaluate  and  analyze  the  performance  of  deriv¬ 
ing  ZiM  from  ocean-color  remote  sensing,  Figure  2  (2a  for 
rPAR  =  0.5  (50%),  2b  for  rf>AH  =  10%,  and  2c  for  rPAR  =  1%) 
and  Figure  3  present  /fra-derivcd  values  (zf^)  versus  those 
determined  from  rPAR{z)  (zj11™),  and  Table  1  summarizes 


their  errors.  Quantitatively,  an  averaged  percentage  error  for 
all  observations  is  calculated  from 


£  = 


(13) 


[21]  Measurements  at  one  station  were  excluded  from 
these  analyses  because  Z|%  from  remote-sensing  reflectance 
(no  matter  which  method)  is  about  four  limes  the  zt%  from 
PAR  profile.  This  one  station  is  regarded  as  an  outlier 
because  of  the  significant  difference  between  in-water  and 
above-water  detenmi nations. 

4.L  IOP-C entered  Approach 

[22]  For  this  data  set  (64  stations)  that  covers  both 
oceanic  and  coastal  waters,  the  average  percentage  error 
(e)  between  modeled  and  measured  property  is  18.5%  (with 
a  maximum  error  of  70.9%)  for  z™%  (for  a  range  of  0,6- 
18.4  in  from  PAR  measurements);  while  the  emors  are 
13.8%  (maximum  error  of  61.2%)  for  zJ0%  (2.1-47.1  m), 
and  13.7%  (maximum  error  of  63.5%)  for  zt%  (4.3  to 
82.0  m).  The  root-mcan-squarc  error  in  log  scale  (RMSE) 
is  0.079  for  z1%1  significantly  smaller  than  the  RMSE 
(0.329)  of  ^-derived  [Chi]  (also  see  Figure  4). 

[23]  For  z$o%f  excluding  the  residual  errors  in  the  Kv  1S 
model  and  the  QAA  algorithm,  there  are  a  few  extra  sources 
to  contribute  to  its  discrepancies.  These  extra  sources 
include  (I)  R4R(z)  attenuates  sharply  at  surface  [Pauimn 
and  Simpson,  1977;  Zaneveld  and  Spinrad ,  1980],  making 
it  harder  to  precisely  determine  z5u*y.  from  the  profile  of 
rPAk{z)\  (2)  surface  layer  suffers  more  influence  from  the 
wavy  surface  [Zaneveld  et  at 2001],  which  can  converge  or 
diverge  incoming  solar  radiation;  and  (3)  variable  ship- 
perturbation  to  the  near-surface  sensor  [Gordon,  1985]. 

[24]  Additionally,  larger  error  came  from  the  two  dear- 
water  stations  (the  two  with  circles  in  Figure  2),  where 
surface  chlorophyll-a  values  were  around  0.07  mg/m3.  The 
two  stations  had  similar  [Chi]  values  and  similar  Rn  spectra, 
and  nearly  identical  zil)%  and  z1%  depths  from  PAR(z) 
profiles  (see  Figures  2b  and  2c);  their  z50%  depths  from 
PAR(z ),  however,  differed  by  nearly  80%.  It  Is  not  clear  yet 
what  might  have  caused  such  big  differences  in  the  z50% 
from  PAR(z)  measurements.  Excluding  these  two  stations, 
the  average  error  for  z50%  dropped  to  16.8%. 

[25]  Much  better  results  are  achieved  for  zlu%  (and  Z|%), 
though,  again,  larger  errors  occurred  at  the  two  clearer- water 
stations  (Figure  2b,  in  circle).  Unlike  z50%,  however,  the 
^10%  values  from  PAR(z)  are  about  the  same  for  these  two 
stations.  Overall,  it  appears  that  there  is  a  trend  of  slight 
underestimation  for  zJ0%  in  the  range  of  10-30  m  for  this 
historical  data  set.  High-quality  measurements  of  both  lOPs 
and  PAR(z)  profiles  are  required  to  isolate  the  error  sources. 

[26]  For  the  entire  range  of  zl%  (4.3-82.0  m)  measured  at 
different  times  and  from  different  regions,  the  RTS  derived 
values  clearly  matched  the  /H/?(z)-delermincd  values  very 
well.  Though  many  potential  sources  of  error  prevent  exact 
agreement  between  the  two  independent  data  sets,  the  low 
average  error  (13.7%  in  linear  scale,  0.079  in  RMSE)  for 
zl%  indirectly  validated  the  approach  of  deriving  lOPs  from 
R re  and  calculating  £vls(z)  from  lOPs.  It  is  even  more 
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Figure  2.  Comparison  between  from  lOP-centered  approach  and  z,%  from  in  situ  PAR  profile  (a)  for 
^50%»  (b)  for  Zio%,  and  (c)  for  Z|%.  The  two  circled  stations  (Figure  2a)  had  nearly  identical  zJ(>% 
(Figure  2b)  and  z1%  (Figure  2c)  values. 


encouraging  given  that  the  measurements  were  not  made  in 
vertically  homogeneous  waters*  If  we  use  the  vertical 
variation  (within  the  euphotic  zone)  of  the  diffuse  attenua¬ 
tion  coefficient  [Kirk,  1994]  at  440  nm  (^d(44G))  as  a  proxy 


to  represent  the  non -homogeneity  of  optical  properties  in 
the  upper  water  column*  the  average  coefficient  of  variation 
is  17%  (maximum  is  51%)  for  this  data  set  If  the  variation 
is  determined  by  chlorophyll  fluorescence  profile,  the  aver- 
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Figure  3.  Comparison  between  from  [Chl]-centcred  approach  and  zl%  from  in  situ  PAR  profile.  z!% 
is  calculated  by  equation  (10),  and  all  waters  are  assumed  Case-]. 


age  coefficient  of  variation  is  23%  (maximum  is  91%)* 
These  values  indicate  varying  optical  or  biological  proper¬ 
ties  in  the  upper  water  column  foF  each  station.  The 
excellent  results  of  remotely  estimated  z1%  and  Z|0%  appar¬ 
ently  indicate  that  lOPs  from  the  semianalytical  algorithm 
and  the  analytical  KwS  model  integrate  and  largely  com¬ 
pensate  for  some  of  the  vertical  variations.  At  least  in  part, 
this  is  because  remote-sensing  reflectance  is  always  a  mea¬ 
sure  of  the  weighted  average  of  optical  properties  in  the  upper 
water  column  [Gordon  and  Clark ,  1980;  Sathyendranath  and 
Platt ,  1989;  Zone  veld  ei  al%  2005]. 

[27]  It  is  important  lo  recognize  that  what  was  carried  out 
here  were  not  regressions  from  the  data  set  to  develop  an 
empirical  relationship,  but  rather  were  independent  tests  of  a 
system  developed  earlier  from  other  sources.  Such  efforts, 
together  with  the  resultant  small  error,  provide  us  much 
higher  confidence  in  the  application  of  such  semianalyti- 
cally  based  approaches  to  broader  ranges  of  environments. 
Additionally,  there  Is  no  requirement  of  the  waters  to  be 
Case-1  for  the  application  of  this  lOP-centered  approach, 
thereby  avoiding  the  hurdle  of  identifying  a  water  body  as 
Case-1  or  not  [Lee  and  Hut  2006;  Mobley  ei  ah,  2004] 
before  processing  the  data.  The  advantage  of  this  IOP- 
c entered  approach  is  further  supported  by  comparing  the 
2 1%  values  that  are  calculated  from  /^-derived  chlorophyll 
concentrations  (see  Figure  3  and  below), 

4*2,  [Chi [-Centered  Approach 

[2s]  Figure  3  shows  z | %  derived  from  /?rs  compared  with 
measured  z  1%  from  PAR{z)  for  the  same  data  set,  but  this 
time  the  concentration  of  chlorophyll-a  is  derived  first 
empirically  from  Rm  and  with  the  assumption  that  the 
waters  under  study  fits  the  Case-1  category  (equations  (10) 
and  (11)).  The  averaged  error  for  zJ%  using  this  approach 
is  32,7%  (maximum  error  is  218%),  whereas  the  average 
error  of  OC4v4  derived  surface  chlorophyll  (Figure  4)  is 
75.2%  with  a  maximum  of  565%.  Generally,  the  derived 
zl%  from  /?rs-[Chl]  matched  the  measured  values  better 


for  waters  with  z1%  deeper  than  ^30  m  (see  Figure  3), 
but  overestimated  zt%  for  depths  shallower  titan  ^30  m. 
This  might  he  because  waters  with  deeper  zJ%  are 
clearer,  providing  a  better  fit  lo  the  Case- 1  assumption, 
although  oceanic  waters  are  not  necessarily  Case-1  [Lee 
and  Hu,  2006;  Mobley  ei  al.t  2004],  Waters  with 
shallower  z1%  were  often  in  coastal  regions,  and  their 
optical  properties  were  less  likely  to  co-vary  with 
chlorophyll  concentration.  This  is  also  consistent  with 
an  earlier  study  regarding  the  down  welling  diffuse  atten¬ 
uation  coefficient  [Lee  ei  ait  2005a],  where  its  value 
was  underestimated  (also  based  on  /?,s-derived  [Chi])  for 
higher  values  (shallower  Z]%),  but  was  quite  good  for 
lower  values  (deeper  zl%).  Compared  with  the  z]% 
estimated  from  OC4v4-[Chl],  estimated  from  the 
JOP-centered  approach  are  generally  shallower  and  more 
consistent  with  measurements,  which  would  lead  to  a 
smaller  compensation  irradianee  if  derived  as  in  Siegel  et 
at  [2002]. 

[29]  It  is  necessary  lo  emphasize  that  "the  zl%  relation¬ 
ship  developed  in  Morel  [1988]  and  Morel  and  Mari  town  a 
[2001]  requires  either  the  mean  chlorophyll  concentration  - 
or  the  water-column  -  Integra  ted  concentration  -  within  the 
euphoric  zone  as  input”  ($.  Maritorena,  personal  commu¬ 
nication),  Because  of  the  existence  of  subsurface  maxima  of 
chlorophyll  concentration,  the  mean  value  is  normally 
greater  than  the  surface  value.  Consequently,  if  surface 
chlorophyll  (e.g.,  the  product  from  current  ocean -color 


Table  I.  Properties  and  Their  Error  Characters  When  Derived 
From  Ocean- Co  I  or  Remote  Sensing 


JGF-Ccntered  Approach 

Chi-Centered 

Approach 

*SMt 

Z|% 

*1% 

fChlj 

£  (average  rrror) 

18,5% 

13,8  V. 

13.7% 

32.7% 

75.2% 

Max.  error 

70,9% 

61.2  V. 

63.5% 

218% 

565% 

RMSE 

0.117 

(1.077 

0.079 

0.162 

0.329 
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Figure  4.  Comparison  between  OC4v4  derived  [Chi]  and  measured  [Chi]  from  surface  samples. 


remote-sensing  algorithms)  is  used  for  the  calculation  of 
it  is  very  likely  that  significant  overestimation  of 
values  will  result,  as  shown  by  Figure  5.  Here  z}%  is 
calculated  simply  using  the  relationship  developed  in  Morel 
and  Maritorena  [2001].  Since  we  do  not  always  know  the 
details  of  the  vertical  distribution  of  chlorophyll  concentra¬ 
tion  for  each  station  (especially  from  remole  sensing),  it  is 
assumed  arbitrarily  dial  the  mean  concentration  of  chloro¬ 
phyll  within  the  euphoric  zone  is  1.3  limes  that  derived  by 
the  OC4v4  algorithm.  With  this  consideration,  the  average 
error  Tor  "i%  is  343%.  The  average  error  is  much  larger  (e  = 
49.9%)  when  no  such  adjustment  is  considered,  but  can  be 


reduced  to  223%  if  the  mean  concentration  is  assumed  as 
1.8  times  the  GC4v4  derived  surface  value. 

5.  Conclusions  and  Implication  to  Ocean-Color 
Remote  Sensing 

[an]  In  this  study,  euphoric  zone  depths  (and  two  other 
environmental-optics  depths)  were  derived  semianalytically 
from  spectral  remote- sen  sing  reflectance  for  measurements 
made  in  the  Arabian  Sea,  the  Gulf  of  Mexico*  and  the 
Monterey  Bay.  The  derived  values  were  then  compared  with 
the  values  from  in- water  vertical  profiles  of  PAR.  ft  was 


0  20  40  60  50 

zl%  from  PAR(z)  [mj 

Figure  5.  Comparison  between  zi%  from  [Chl]-centered  approach  and  z1%  from  in  situ  PAR  profile* 
Here  zl%  is  calculated  by  rearranging  equation  (6)  of  Morel  and  Maritorena  [2001],  and  the  average  [Chi] 
in  the  euphoric  zone  is  assumed  as  13  times  the  surface  value  derived  by  OC4v4  (see  text  for  details)* 
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found  that  the  two  agreed  with  each  other  very  well  for  the 
entire  range  of  the  data  set.  The  average  percentage  error  in 
linear  scale  is  18,5%  for  z50%  (in  a  range  of  0.6-18.0  m), 
13.8%  for  zl0%  (2.1-47,1  m),  and  13.7%  for  z!%  (4, 3“ 
82.0  m).  Such  small  errors  suggest  a  closure  between  Hie 
two  independent  measurements  and  determinations,  and 
indirectly  validate  the  semianalytical  derivation  of  lOPs 
from  /frs  and  the  lOP-centered  model  of  AVis-  To  ensure 
its  reliable  applications  to  broad  range  of  waters,  however, 
more  tests  and  validations  with  a  wider  dynamic  range  are 
certainly  desired,  especially  for  waters  with  zt%  deeper  than 
100  m. 

[it]  Also  presented  are  comparisons  of  z1%  calculated 
from  the  same  R but  based  on  a  methodology  that  uses 
chlorophyll  values  and  the  Case-1  water  assumption.  The 
average  error  is  much  larger  (32.7%)  when  compared  with 
Zl%  from  measurements,  with  better  results  for  oceanic 
waters,  but  overestimates  Zi%  for  coastal  waters.  As  dis¬ 
cussed  earlier,  the  primary  source  contributing  to  this  larger 
error  includes  the  inclusion  of  coastal  waters  with  many 
stations  hardly  belonging  to  the  Case-1  category.  An  over¬ 
age  error  of  75%  in  [Chi]  determined  by  the  OC4v4 
algorithm  for  these  waters  suggests  that  use  of  a  global 
algorithm  for  coastal  waters  imparts  a  significant  contribu¬ 
tion  to  the  error  field, 

[32]  The  lOP-centered  approach,  however,  worked  not 
only  with  oceanic  waters,  but  also  with  the  more  complex 
coastal  and  shelf  waters,  as  represented  by  the  data  collected 
in  this  study.  Also,  because  the  relationships  between  zx% 
and  Rrs  are  linked  analytically  by  IOPs,  this  IGP-centered 
approach  avoids  parameterization  regarding  the  wide  re¬ 
gional  and  seasonal  bio -optical  variations.  This  is  extremely 
important  if  an  algorithm  is  going  to  be  applied  to  get 
quantitative  global  observations.  In  the  current  process  of 
deriving  [Chi]  of  global  oceans  from  satellite  ocean-color 
remote  sensing,  the  “  global  “  algorithm  is  with  a  single  set 
of  parameters  that  does  not  correct  for  regional  and  seasonal 
bio-optical  variations  when  applied  lo  global  oceans.  If 
proper  regional  and  seasonal  p ammeter izati on s  are  utilized 
as  desired  [Carder  el  al,  1989:  O'Reilly  et  a l .,  1998],  a 
different  global  mean  value  for  chlorophyll,  and,  likely,  a 
different  temporal  trend  of  this  mean  relative  to  those  of 
Gregg  and  Conkrighi  [2002]  and  Antoine  et  at.  [2005], 
would  be  reached  (Shang  and  Lee,  manuscript  in  prepara¬ 
tion,  2007). 

[b]  EnphoSic  zone  depth  (z{%)  (or  1 0%-light  depth,  zJ0%) 
represents  depths  where  only  1%  (or  10%)  of  surface  PAR 
remains.  Compared  with  the  measurement  of  chlorophyll 
concentration  (large  uncertainties  exist  between  measured 
by  FIPLC  method  and  measured  by  fluorometric  method 
[Muelter  et  at.,  2002b;  Trees  et  at.,  1985]),  or  measurement 
of  absorption  and  scattering  coefficients  of  the  bulk  water  or 
phytoplankton,  z1%  (or  zl0%)  is  much  easier  and  more 
accurate  to  determine  in  the  field.  It  does  not  require  an 
absolute  radiometric  calibration  of  the  PAR  sensor,  and  it 
does  not  require  delicate  calibration  of  the  final  product  or 
complicated  and  commonly  error-bearing  post  processing 
(e.g.,  the  absorption  coefficient  of  phytoplankton  from 
filter-pad  technique  [Aitati  et  at .,  1995;  Cleveland  and 
f Veidemann,  1993;  Mitchell,  1990]). 

[34]  Also,  the  close  agreement  between  tfrs -calculated  and 
field-measured  z1%  (and  z10-/J  can,  in  large  part,  be  attrib¬ 


uted  to  the  fact  that  both  are  cumulative  measurements  of 
the  upper  water  column.  Chlorophyll  concentrations  or 
absorption  and  back  sea  tiering  coefficients  from  discrete 
water  samples  may  not  accurately  represent  the  average  of 
the  upper  water  column  unless  it  is  uniform;  and,  conse¬ 
quently,  larger  uncertainty  will  be  introduced  when  com¬ 
paring  the  values  from  discrete  water  samples  with  those 
from  remote  sensing.  Many  of  those  uncertainties  would 
be  minimized  for  z[%  (or  ri0«/B),  and  we  can  safely  set  an 
accuracy  goal  of  within  20%  from  satellite  measurements 
after  further  evaluation  and  refinement.  Remote-sensing 
products  with  such  a  small  error  will  bring  confidence  on 
quantitative  ocean-color  remote  sensing  and  boost  inter¬ 
ests  and  applications  of  professionals  and  the  general 
public. 

[35]  Because  of  these  distinctive  and  unique  character¬ 
istics,  and  because  the  error  of  measured  z1%  is  significantly 
smaller  than  that  of  chlorophyll-a  concentration,  zl%  (and/or 
*  1 0%)  can  be  an  ideal  parameter  to  validate  a  system  of 
ocean-color  remote  sensing.  More  importantly,  z1%  (or  zl0%) 
measures  water  clarity  much  more  rigorously  than  Secchi 
depth  [ Preisendor/er ,  1986]  and  provides  much  more  rcli* 
able  results.  The  variations  of  zt %  (or  zy>%)  (after  removing 
sun  angle  effects)  indicate  clearly  changes  of  water  quality 
of  an  ecosystem.  A  time  series  ofzl%  (or  r10%)  of  the  global 
oceans  from  satellite  measurements,  combined  with  histor¬ 
ical  and  new  measurements,  will  provide  us  unprecedented 
and  confident  evaluation  of  the  distribution  and  trend  of 
water  clarity  of  the  world  oceans.  Note  that  water  clarity  has 
profound  effects  on  primary  production  [Behrenfeld  and 
Falkowski,  1997;  Platt  and  Sathyendranath,  1988]  and  heat 
deposition  [Kirk,  1988;  Lewis  et  at 1990;  Morel  and 
Antoine,  1994;  Zone  veld  et  ai ,  1981]  in  the  upper  water 
column. 

Appendix  A:  Updated  Quasi-Analytical 
Algorithm  (QAA_v4) 

[36]  The  Quasi- Analytical  Algorithm  (QAA)  was  devel¬ 
oped  by  Lee  et  at.  [2002]  to  derive  the  absorption  and 
buckscattering  coefficients  by  inverting  spectral  remote¬ 
sensing  reflectance  (/?,*(  A)).  QAA  starts  with  the  calculation 
of  the  total  absorption  coefficient  (a)  at  a  reference  wave¬ 
length  (Ao)>  and  then  propagate  the  calculation  to  oilier 
wavelengths.  To  briefly  summarize,  this  algorithm  consists 
of  the  following  elements: 

[3?]  1 .  The  ratio  of  backscattcring  coefficient  (bh)  to  the 
sum  of  backscattcring  and  absorption  coefficients  + 
frb))  at  A  is  calculated  algebraically  based  on  the  models  of 
Gordon  et  at.  [1988]  and  Lee  et  at.  [1999], 

M  A)  —Q.OS95  +  /0.008  +  0A99r^{\) 
a(\)  +  bh(\)  0.249  fA  * 

Here  rr*(A)  is  the  spectral  remote-sensing  reflectance  just 
below  the  surface  and  is  calculated  from  /?rv(A)  through. 


r«(  A)  =  *,(  AX0.52  +  \JR„(\))-  (A2) 
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Figure  Al.  Comparison  between  IGPs  derived  from  R ^  and  lOPs  from  in  situ  measurements.  QAA_v4 
was  used  for  the  derivation  of  lOPs,  and  Rzs  values  at  410,  440,  490,  555,  and  670  mil  were  used  as 
inputs. 


[as]  2<  The  spectral  hb(X)  is  modeled  with  the  widely  used 
expression  [Gordon  and  Morel,  1983;  Smith  and  Baker, 
1981], 

=  +  (A3) 

where  hbw  and  hb p  are  the  baekscattering  coefficients  of 
pure  seawater  and  suspended  particles,  respectively.  Values 
of  bbyf( A)  are  provided  in  Morel  Jl 974]. 

[37]  3.  When  n(A0),  the  ratio  of  bj{a  +  bh)  at  \h  and 
^bw(\)  arc  known,  ^(Aq)  in  equation  (A3)  can  be  easily 
derived  with  the  combination  of  equations  (Al)  and  (A3). 
The  values  of  6bp( A)  at  other  wavelengths  are  then  calcu¬ 
lated  after  the  power  parameter  (?/)  is  estimated  from  [Lee  et 
aiy  2002], 

(M) 

[4n]  4.  Finally,  by  applying  /?b( A)  to  the  ratio  of  b^J(a  + 
hb)  at  A  (equation  (Al)),  the  total  absorption  coefficient  at  A, 
n(A),  is  calculated  algebraically. 


[41]  The  update  of  die  QAA  is  related  to  the  calculation 
of  n(Ao).  In  this  updated  version  of  QAA  (v.4)t  it  has  two 
selections  for  the  calculation  of  a(555):  one  for  data  sets 
where  there  are  no  Rn  measurements  at  longer  wavelengths 
(640  or  670  nm,  for  instance);  and  one  for  data  sets  where 
there  are  measurements  at  those  longer  wavelengths. 
And  the  second  round  of  calculation  [Lee  et  al,  2002]  is 
removed. 

[42]  For  data  sets  having  no  at  longer  wave¬ 
lengths,  a(555)  is  now  estimated  using  the  ratio  p  = 
log(max(/?,.7(44Q,  490,  51 0)Vi?,v(555))  (for  use  with  the 
SeaWiFS  sensor),  and  takes  tire  form, 

K(555)  -  0.0605  +  [Q~l -i3V+o.4iV-fl*«V(  (A5) 


[43]  The  constants  in  equations  (A5)  and  (A6)  arc  combi¬ 
nations  of  those  from  OC4v4  [O'Reilly  et  al,  2000]  and 
Morel  and  Maritorena  [2001], 

[*14]  For  data  sets  having  R r*  at  640  nm  (note  that  R r3t  at 
longer  wavelengths  are  very  important  for  inversion  of 
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waters  wilJi  high  absorption  coefficients  [ IOCCG ,  2006;  Lee 
el  al.,  2005a]),  o{555)  is  now  estimated  as  follows. 


X  =  I°E 


/  *,-,(440)  +  /;,, (490)  \ 

^„(555)  +  2««„(640)J’ 


(A7) 


a(555)  =  a,„{555)  +  (AS) 

with  constants  in  equation  (AS)  obtained  by  least  squares 
fitting  tf(555)  of  the  synthetic  data  set  adopted  by  the 
IOCCG  [IOCCG,  2006]. 

[45]  In  the  original  version  of  QAA  [Lee  et  ait  2002]  a 
relatively  arbitrary  boundary  was  used  for  the  transition 
between  values  calculated  with  555  nm  and  640  nm  as 
reference  wavelengths,  respectively.  In  this  updated  QAA 
(equations  (A7)  and  (A8)),  now  seamless  transition  is 
achieved  for  waters  from  open  ocean  to  coastal  region. 
Available  measurements  suggest  that  QAA_v4  is  applicable 
to  waters  with  #(440)  as  high  as  ^2.0  m~K  Simply  change 
the  corresponding  values  at  555  nm  to  that  at  550  nm, 
QAA_v4  (equations  (A7)  and  (AS))  can  be  applied  to  data 
from  both  SeaWiFS  and  MOD  IS  without  significant  effects 
to  the  final  results. 

[4^]  For  tlie  SeaWiFS  and  MODiS  sensor,  there  is  no 
band  to  measure  R ^  at  640  nm,  but  there  are  bands  to 
measure  /?rs  at  667  and/or  678  nm,  For  such  cases,  a 
simulated  /?ts(640)  needs  to  be  generated  as  in  Lee  et  ai 
[2005a],  which  is 

40)  -  0.01  Rr](555)  +  l-4/J„(667) 

-  0.0005^(667)7^(490),  (A9) 

and  this  simulated  /?rs(640)  is  kept  greater  than  L2  R^( 667), 

[4?]  Figure  Al  compares  optical  properties  derived  from 
RIS  by  QAA_v4  (a(555)  from  equation  (A8))  with  that  from 
water-sample  measurements,  where  data  is  from  the 
NOMAD^OCBAM  data  set  [WerdeU  ttnd  Bailey*  2005], 
but  subsampled  for  valid  measurements  of  both  /?rs  and 
lOPs  (T.  Smyth,  personal  communication).  Clearly,  excel¬ 
lent  agreements  were  achieved  between  the  two  indepen¬ 
dent  determinations.  To  indicate  the  applicability  of 
QAA_v4  with  satellite  dala,  it  was  the  simulated  /?rs(640) 
(equation  (A9))  instead  of  the  measured  /frs(640)  used  in  the 
derivation  of  lOPs  (and  the  zt%)  in  this  study. 

[48  ]  For  a  recent  measurement  made  in  the  south  Pacific 
gyre  (M.  Lewis,  personal  communication),  consistent 
results  (with  zl%  around  150  m)  were  obtained  from 
measurements  of  the  vertical  profiles  and  the  remote  sensing 
reflectance. 
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